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Abstract. We study the relationship between the multiplicity of 
a fixed point of a function <?, and the dependence on e of the length 
of e- neighborhood of any orbit of g, tending to the hxed point. The 
relationship between these two notions was discovered in Elezovic, 
Zubrinic, Zupanovic 5J in the differentiable case, and related to 
the box dimension of the orbit. 

Here, we generalize these results to non-differentiable cases in- 
troducing a new notion of critical Minkowski order. We study 
the space of functions having a development in a Chebyshev scale 
and use multiplicity with respect to this space of functions. With 
the new definition, we recover the relationship between multiplic- 
ity of fixed points and the dependence on e of the length of e- 
neighborhoods of orbits in non-differentiable cases. 

Applications include in particular Poincare maps near homo- 
clinic loops and hyperbolic 2-cycles, and Abelian integrals. This is 
a new approach to estimate the cyclicity, by computing the length 
of the e-neighborhood of one orbit of the Poincare map (for exam- 
ple numerically), and by comparing it to the appropriate scale. 
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1. Introduction 

The multiplicity of a fixed point of a differentiable function can be 
seen from the density of its orbit near the fixed point as was shown 
in [5] . We recall this result in Theorem [I] The information on the 
density is contained in the behavior of e— neighborhood of the orbit 
near the fixed point and is usually measured by the box dimension of 
the orbit. It was further noted in [HI] that the box dimension of the 
orbit of Poincare map around a focus or a limit cycle shows how many 
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limit cycles can appear in bifurcations. This gave an application of the 
result from Theorem [I] to continuous dynamical systems. 

The idea of this article is to generalize these results to a class of 
functions which are non-differentiable at a fixed point. The goal is 
again to estimate the multiplicity of a fixed point of such a function 
only from the asymptotic behavior of the length of the e-neighborhood 
of any of its orbits close to the fixed point, as e — >■ 0. The results 
can be applied to continuous dynamical systems. While differentiable 
functions described above appear as displacement functions near limit 
cycles and foci, non-differentiable functions appear naturally as dis- 
placement functions near polycycles, see e.g. [H], [15] (see Section [1]). 
The multiplicity of a fixed point of the displacement function near 
some limit periodic set reveals the number of limit cycles that appear 
in the unfoldings of the limit periodic set. It is of interest to find at 
least an upper bound on the multiplicity. 

Calculating numerically an orbit of the Poincare map of the limit 
periodic set, one can estimate the length of its e-neighborhood for small 
values of e and thus estimate its asymptotic behavior. 

In the differentiable case (foci, limit cycles), see TheoremjlJ it suffices 
to compare the behavior of the length with discrete scale of powers, 
{e, e 1//2 , e 1//3 , . . .}. The moment when comparability occurs reveals 
multiplicity, i.e. cyclicity. We see additionally that this moment is 
signaled by the limit capacity (the box dimension) of the orbit which 
actually shows the density of the orbit around the fixed point: the 
bigger this density is, more limit cycles can appear in perturbations. 

In non-differentiable cases, however, we show that it is not sufficient 
to compare the length of the e-neighborhood with the scale of powers to 
estimate multiplicity and cyclicity. The idea behind Theorems [2] and [3] 
is in finding the appropriate scale to which the length should be com- 
pared to obtain precise information on the multiplicity. This scale, 
as we will see, depends on the unfolding and should be estimated at 
least from above. Here, instead of box dimension, the new notion of 
critical Minkowski order is introduced, to signal the moment when the 
comparability occurs in the new scale. 

The article is organized as follows. First, in Subsection 1.1 we recall 
the connection from [5] between the box dimension of the orbit and the 
multiplicity of the fixed point in the differentiable case, see Theorem [T] 
In Subsection 1.2 we recall and introduce definitions and notions we 
need in non-differentiable cases. Finally, in Section [2j we state our 
main results concerning non-differentiable cases, see Theorem [2] and 
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Theorem [3j Some applications to continuous dynamical systems are 
given in Section |4} 

1.1. Differentiable case. 

Denote Diff r [0,d) the space of C r -differentiable functions on [0,d), 
for r sufficiently big, d > 0. Let / G Diff r [0,cf), /(0) = and 
x > f(x) > 0, for x G (0, d). Put 

(1) g = id-f 

and consider the orbit S 9 (xo) of < xq < d by g: 

(2) ^(xo) = {x n | n G N}, x n+ i = g(x n ). 

Let jUg (g) be £/ie multiplicity ofO as a fixed point of the function g 
in the family Diff r [0, d). That is, the number of fixed points that can 
bifurcate from by bifurcations within Diff r [0, d). Then, 
(3) 

Vi*(g) = k, if /(0) = /' (0) = . . . = /(^(O) = 0, /«(0) ^ 0, 

i.e., is a zero of multiplicity /xo(/) = k of /. 

Now we define the Minkowski content and the box dimension of a 
bounded set. Let U C M. N be a bounded set. Denote by |A e (C/)| the 
Lebesgue measure of e— neighborhood of U. 

By lower and upper s- dimensional Minkowski content of U, s > 0, 
we mean 



M'JU) = liminf \ A i U ^ and M* S (U) = limsup 

respectively. Furthermore, lower and upper box dimension of U are 
defined by 



&m p U = inf{s > | Mt(U) = 0}, dim B U = inf{s > | M* S (U) = 0}. 

As functions of s G [0,N], M* S (U) and Ml{U) are step functions that 
jump only once from +oo to zero as s grows, and upper and lower box 
dimension contain information on jump in upper and lower content 
respectively. 

If d\m B U = dirngt/, then we put dims([/) = dhn B U = dim^C/ and 
call it the box dimension of U . In the literature, upper box dimension 
(also called limit capacity) has been widely used. For more details on 
box dimension, see Falconer |6| or Tricot |16|. 



Our case is 1-dimensional so iV = 1 in the definition of box dimen- 
sion, and also U C [0, d), where d > 0. We are interested in measuring 
the density of accumulation of the orbit of a function g near its fixed 
point zero. Let g be sufficiently differentiable on [0,d), d > 0, such 
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that g(0) = 0, we denote by S 9 (xq), < xq < d, the orbit of xq by g 
defined by x n+ \ = g(x n ), xq < d, and tending monotonously to zero. 
In 1-dimensional differentiable case it is verified that dimes' 9 (xo)) is 
independent of the choice of the point xq in the basin of 0. Therefore 
one can define box dimension of a function g by 

dining) = dim B (S 9 (x )), 

for any xq from the basin of attraction of 0. 

For two positive functions F(x) and G(x), with no accumulation 
of zeros at x = 0, we write F(x) ~ G(x), as x — > 0, if there exist 
two positive constants A and B and a constant d > such that A < 
F(x)/G(x) < B, x G (0, d), and call such functions comparable. In the 
sequel, we write F{x) = o(x), if lim^o = 0. 

Now we reformulate Theorem 1 from [5], connecting box dimension 
and multiplicity in the differentiable case. 

Theorem 1. Let f be sufficiently differentiable on [0,d) and positive 
and strictly increasing on (0,d). Let g —id—f and suppose that the 
multiplicity of as a fixed point of g is finite and greater than 1. That 
is, 1 < /io' x (g) < oo. Let x £ (0, d), S 9 {xq) be defined as in ^ and 
let \A e (S 9 (xo))\ be the length of the e -neighborhood of the orbit S 9 (xq), 

£>0. 

Then 

(4) \A £ {S 9 {x ))\ -s 1 ' ^ X(9 \ as e^0. 

If fiQ lx (g) = 1 and additionaly f(x) < x on (0, d), then 

(5) \A(S 9 (x n ))\~ { £ (- lo S £ )' i ff'(°)< 1 as£ ^Q 
Moreover, for 1 < ^Q X {g) < oo, 

(6) ^ = l-di!n B (,)- 

Sketch of proof. We illustrate the proof on the simplest case when g is 
linear, g(x) = kx, k G (0,1). Take any initial point x G (0, d). By 
recursion, it is easy to compute the whole orbit by g: 

(7) x n = k n x , neN. 

To compute the asymptotic behavior of the length of the 
e-neighborhood of the orbit, we divide the e- neighborhood in two parts: 
the nucleus, N e , and the tail, T £ . The tail is the union of all disjoint 
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(2£)-intervals of the e-neighborhood, before they start to overlap. It 
holds that 

(8) \A £ (S°(x ))\ = \N e \ + \T £ \. 

Let n £ denote the index separating the tail and the nucleus. It de- 
scribes the moment when (2e)-intervals around the points start to over- 
lap. We have that 

(9) \N £ \=x nE +e, |T e | ~n £ -e, e 0. 

To find the asymptotics of n £ , we have to solve x Ue+ i — x ne ~ 2e, e — > 0. 
Using Q, we get 

x ne ~ £, n £ ~ - log £, e ->• 0. 

From (|9]), we get 

|JV e | ~ e, \T £ \ ~ e(-loge), 
therefore, by g, \A £ (S g (x ))\ ~ e(-loge), as e -> 0. □ 

Note that in Theorem [T] we assume / to be differentiable at zero 
point x — 0. In this article, we generalize Theorem [T] to some non- 
differentiable functions at x — 0. 

Since in the non-differentiable cases standard multiplicity of zero is 
not well defined, we use the notion of multiplicity of a point as zero of 
f(x) with respect to a family of functions (see Definition [TJ. The family 
of functions we consider will be the family of functions having a finite 
codimension asymptotic development with respect to a Chebyshev scale 
(see Definition [2]). 

As the main results, in Section [2] we extend formula ^ to non- 
differentiable case. Therefore we have to introduce the notion of criti- 
cal Minkowski order with respect to a Chebyshev scale, see Definition [5j 
This notion is in the differentiable case directly related to box dimen- 
sion, see Remark |2}m). In non-differentiable cases, it generalizes the 
notion of box dimension in a way that a formula similar to formula ^ 
holds. 

1.2. Non-differentiable cases. 

Let us recall some definitions we use in the non-differentiable cases. 

Definition 1. Let A be a topological space and let T = {f\\ A G A} ; 
fx : [0, d) — > E ; be a family of functions. Let Ao G A ; we say that x = 
is a zero of multiplicity greater than or equal to m of the function f\ 
in the family of functions J 7 if there exists a sequence of parameters 
A n — > A , as n — >■ 00 , such that, for every n G N, f\ n has m distinct 
zeros y™, . . . , y^ G [0, d) different from x = and y™ — > ; as n —¥ 00, 
j = l,...,m. 
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We say that x = is a zero of of multiplicity m of the function f\ 
in the family T and write fioifxojJ 7 ) = m, if m is the biggest possible 
integer such that the former holds. 

Putting gx = id — f\, Q = {g\ = id — fx}, the multiplicity of 
as a fixed point of gx with respect to the family Q is iil lx (gx , Q) = 

In [T5] Roussarie introduced the notion of cyclicity, measuring the 
number of limit cycles (isolated periodic orbits) that can be born from 
a certain set called limit periodic set by deformation of a given vector 
field in a family of vector fields. In the case of the family of Poincare 
maps of the family of vector fields in a neighborhood of a limit periodic 
set, cyclicity is given by the above notion of multiplicity of a fixed point 
zero of gx in the family of Poincare maps. Due to the possible loss 
of differentiability of the family of Poincare maps near a limit periodic 
set, the more general notion of multiplicity introduced in Definition [T] 
is necessary. 

Remark 1 (Relation between classical multiplicity and multiplicity in 
a family for differentiable functions). Note that if f is differentiate, 
f G Diff r [0,d), then the classical notion of multiplicity of as a zero of 
f , /io(/) < r as in measures the maximal number of zeros of fx that 
can appear near 0, for fx G T close to f , where T = Diff r [0, d) and the 
distance function is given by d(f,g) = sup fe=0j r |/^(0) — g^(0)\, i.e. 
A*o(/) — Hoi f i Diff r \fl , d)) . Note that the number of zeros /i (/, J 7 ) that 
can appear by deformations J 7 really depends on the family T . Taking 
a family J 7 of deformations bigger or smaller than Diff r [0, d), it is easy 
to give examples with fMo(f, J 7 ) bigger or smaller than Ho{f)- (See e.g. 
Example 1.1.1 and Example 1.1.2 in [T3]). 

We want to study non-differentiable functions having a special type 
of asymptotic behavior at x — 0. The definition of the following se- 
quence of monomials and its properties is based on the notion of Cheby- 
shev systems, see [T2] and [T3], and the proofs therein. A similar notion 
of asymptotic Chebyshev scale is mentioned in [3]. 

Definition 2. A finite or infinite sequence X = {uq, u%, U2, ■ ■ ■} of func- 
tions of the class C[0, d) H Diff r (0, d), r G NU{oo}, is called a Cheby- 
shev scale if: 

i) A system of differential operators D i; i = 0, . . . ,r, is well de- 
fined inductively by the following division and differentiation 
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algorithm: 

D (u k ) = — , 
Wo 

(Di(u k ))' 
D i+1 (u k ) = —— r-, i G N , 

/or ever?/ I: 6 No, except possibly in x = to which they are 
extended by continuity. 

ii) The functions Di{u i+ i) are strictly increasing on [0,d), i G N . 

iii) \im x ^. DjUi(x) = 0, for j < i, % G No- 

We call Di(f) the i—th generalized derivative of / in the scale X. 

Definition 3. A function f has a development in a Chebyshev scale 
X = {u , . . . , Uk} of order k if 

k 

(10) f(x) = y^ayu^z) +ip(x), oti G K, 

t=0 

and t/ie generalized derivatives Di(ip{x)), i = 0, . . . , k, verify Di(i[)(0)) = 
(in the limit sense). 

Note that (in the limit sense) Di(f)(0) — ctj, i — 0, . . . , k. 

Consider a family J 7 = {f\\X G A} of functions having a uni- 
form development of order k in a family of Chebyshev scales T\ = 
(u {x,X),...,u k (x,X)), i.e. 

fc 

(11) f x (x) = ai(\)ui(x, A) + i/>(x, A), A G W. 

i=0 

The development is uniform in the sense that all generalized derivatives 
Djfx, j = 0, . . . , k, can be extended by continuity to x = uniformly 
with respect to A, and this extension is continuous as function of A. 

The following lemma generalizes Remark [Tj It gives the connection 
between the index of the first nonzero coefficient in development of a 
function / G T in a Chebyshev scale and the multiplicity of as a zero 
of / in the family T '. 

Lemma 1. Let T\ = (uo(x, A), . . . , Uk(x, A)), A G A, be a family of 
Chebyshev scales and T = (f\) a family of functions having a uniform 
development in the family of Chebyshev scales X A of order k and f\ Q G 
J- ' . If the generalized derivatives satisfy 

(12) A(/a o )(0) = 0, i = 0, . . . , k - 1, and D ko (f Xo )(0) ^ 0, k < k, 
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i.e., if afc (Ao) is the first nonzero coefficient in the development of f\ , 
then the multiplicity of as zero of f\ in the family J 7 is at most ko 
'i.e., noifx^J 7 ) < k ) 

X '" ' ^ [ox~)i=0...k -l, j=l...k {^0) 



If moreover A C E JV , k < N, and the matrix (^ L )i=o...k -i, j=i...k (\ 



is of maximal rank (i.e., equal to k ), then (12) is equivalent to 

Proof of Lemma Q] is based on Rolle's theorem and the observation 
that dividing by nonzero functions the number of zeros is unchanged. 
If the matrix (|^) is of maximal rank, then by the implicit function 
theorem one can consider parameters Xj, j = 1, . . . , k , as functions of 
at , ■ ■ ■ , Qffco-i, Afco+i, . . . , Ajv and a , ■ ■ ■ , a ko -i, Xk +i, . . . , Ajv as new pa- 
rameters. Then making a sequence of small deformations starting with 
a!jfc _i, . . . ,ak -2 etc, one can create k Q small zeros in a neighborhood 
of 0. For the details of the proof, see Example 1.1.3 in [TBI . 



Example 1. (Examples of Chebyshev scales on [0,d)) 

i) differentiable case: e.g. X = {1, x, x 2 , x 3 , x 4 , . . .}, 
ii) non- differentiable case: 

- X = {x ao ,x a \x a2 , . . .}, cti e R, < a < ati < a 2 < . . . 

- X = {1, x(— logx),x, x 2 (— logx), x 2 ,x 3 (— logs), X 3 , . . .} 

- More generally, X can be any set of monomials of the type 
x k (—logx) 1 , ordered by increasing flatness: 

x l (— log a;)- 7 < x k (—lnx) 1 if and only if (i < k) or (i = j and j > I). 

in) For more general examples corresponding to Poincare map at a 
homoclinic loop see Section \4-l.S\ 



Definition 4. A function f is weakly comparable to powers, if there 
exist constants m > and M > such that 

(13) m < x ■ (log f)'(x) <M,xe (0, d). 

We call the left-hand side of ( fl~3| ) the lower power condition and the 
right-hand side the upper power condition. A function f is sublinear if 
it satisfies lower power condition and m > 1. 

A similar notion of comparability with power functions in Hardy 
fields appears in literature, see Fliess, Rudolph [7] and Rosenlicht [H]. 
A Hardy field if is a field of real-valued functions of the real variable 
defined on (0, d), d > 0, closed under differentiation and with valuation 
v defined in an ordered Abelian group. Let /, g G H be positive on 
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(0,d) and let \im x ^ f(x) = 0, \im x ^ g(x) = 0. If there exist integers 
M, A^gN and positive constants a, (3 > such that 



(14) 



f(x) < ag(x) M and g{x) < f3f(x 



it is said that / and g belong to the same comparability class/are 
comparable in H. 

Let us state a sufficient condition for comparability from Rosenlicht 
|14j . Proposition 4: 

Proposition 1 (Proposition 4 in [2]). Let H be a Hardy field, f(x), g(x) 
nonzero positive elements of H such that lim^o f( x ) — 0; 



lim. 



0. // 



(15) K0°g/)') = K0°gff)'), 

then f and g are comparable. 



(16) 



The condition (15) is equivalent to (see Theorem in [T4"] ) 
(log /)'(*) 



lim 

x^o [log gy(x) 



L, < L < oo. 



Rosenlicht's condition ( 15 ), i.e. ( 16 ) is stronger than our condition (j 1 3[) 
of weak comparability to powers. If lim^o ^ log 1 ^ ^ 



L, < L < oo 



((13) obviously follows), then / is comparable to power functions in 
the sense (14). 



Note that condition (13) excludes infinitely flat functions, but is 
nevertheless not equivalent to non-flatness. If / is infinitely flat (in the 
sense that all its derivatives tend to zero as x — > 0), then it can easily 
be shown by L'Hospital rule that lim^o = 0, for every a > 
and, as a consequence, the inequality (13) cannot be satisfied. The 
contrary is not true. There exist functions that are not infinitely flat, 



but nevertheless do not satisfy (13). As an example, see the Example |3 
in Appendix. 

Example 2 (Weak comparability to powers and sublinearity) . 

i) Functions of the form 



X 



\ogx)P, a > 0, (3 e 



are weakly comparable to powers. 

This class obviously includes functions of the form x a , x a (— logx) 13 
and rz^^p , for a > and (3 > 0. If additionally a > 1, they are also 
sublinear. 
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ii) Functions of the form 



f(z) = — ■ rg, (3 > 0, 



do not satisfy the lower power condition in (13) 
Hi) Infinitely flat functions of the form 

f(x) = e~^, a > 1, 



do not satisfy the upper power condition in (13), but they are sublinear. 

2. Main results 

We now state the main theorems and their consequences. 

Theorem 2. Let f 6 Diff r (0,d) be continuous on [0,d), positive on 
(0,d) and let /(0) = f'(0) = 0. Assume that f is a sublinear function. 
Put g = id — f and let S 9 (xo) = {x n \ n e N} be an orbit of g, xq < d. 
The following formula for the length of the e -neighborhood of the orbit 
S 9 (xo) holds 

(17) \MS\x ))\^f-\e). 

The sublinearity condition m > 1 in the lower power condition can- 
not be omitted from Theorem [2] For counterexample, see Remark [4] in 
Appendix. 

The following definition is a generalization of box dimension in non- 
differentiable case, according to a given Chebyshev scale. There exists 
in the literature the notion of generalized Minkowski content, see He, 
Lapidus [9], and Zubrinic, Zupanovic |18J . It is suitable in the situation 
where the leading term of |A £ ({7)| does not behave as a power func- 
tion, and we introduce some functions usually called gauge functions. 
Driven by the result of Theorem [2j we follow the idea and define the 
generalized Minkowski content with respect to a family of gauge func- 
tions. By Theorem [2j the e-neighborhood \A £ (S 9 (x )) \ should be com- 
pared to the family obtained by inverting the given Chebyshev scale, 
{u^ 1 (e) , u^ 1 (e) , . . .}. Comparing it to the powers of e as in the stan- 
dard definition of the Minkowski content does not give precise enough 
information. Next we define critical Minkowski order which is close to 
the notion of box dimension. Its purpose is to contain information on 
the jump in the rate of growth of the length of e— neighborhood of an 
orbit. 

The upper (lower) generalized Minkowski content defined in Def- 
inition [5] below can be viewed as function of i, i = !,...,£. By 
Lemma[2j2)(6) in Section [3l it is a discrete function which jumps only 
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once from +00 to through some value < M < +00. This is a behav- 
ior analogous to the behavior of the standard upper (lower) Minkowski 
content as function of s, see Section [l] and e.g. [6]. 

Let X = {uq,u\, . . .} be a Chebyshev scale such that monomials u, 
are positive and strictly increasing on (0, d), for i > 1. Suppose that 
/ has a development in the Chebyshev scale of order I on [0, d) and 
moreover that / satisfies assumptions from Theorem [2] and upper power 
condition. Let g — id — f. We have the following definition: 

Definition 5. By lower (upper) generalized Minkowski content of 
S 9 (xq) with respect to Ui, i = 1, . . . , £, we mean 



M^(S 9 (x ),u,i) = liminf 



\MS"(x ))\ 



M*(S°(x ), Ui ) = limsup^^^ 

e^o ^ (e) 

respectively. It can be easily seen that the behavior of the 
e— neighborhood and thus the definition is independent of the choice 
of the initial point xq < d in the basin of attraction ofO. Therefore we 
define 

m(g,l) = max{i > 1 | M*(S 9 (x ), Ui) > 0}, 
m{g,l) = max{i > 1 | M*(S 9 {x ), m) > 0} 

as the lower (upper) critical Minkowski order of g with respect to the 
scale I, when a jump in lower (upper) generalized Minkowski content 
occurs. If m(g,Z) = m(g,I), we call it critical Minkowski order with 
respect to the scale I and denote m(g,X). 

Remark 2 (Box dimension and critical Minkowski order). 

(i) Using Lemma^ it can easily be seen that the upper and lower 
generalized Minkowski contents M.(S 9 (xo) , Ui) pass from the value 
+00, through a finite value and drop to as i grows. Moreover, 
the critical index i is the same for upper and lower content and 
therefore m(g,I) = iq. 

(ii) If f G Dif f r [0,d) is a differ entiable function, then it has an 
asymptotic development in the differentiable Chebyshev scale, 
X = {1, x, x 2 , . . . x r } of order r. The box dimension and critical 
Minkowski order are then directly related by the formula 



(18) 
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Indeed, assume f(x) ~ x k , 1 < k < r. By Q, \A £ (S 9 (x ))\ ~ 
e 1 ' 1 ". This gives m(g,T) = k. On the other side, by definition, 
the box dimension is the value s such that 1/k — 1 — s. 



(iii) By analogy with the differentiate case (18), we can define gen- 
eralized box dimension of a function g with respect to a Cheby- 
shev scale X by 

i 

dim GB (g,X) = 1 - 



m(g,l) 



This definition is obviously independent of xo from the basin of 
attraction ofO. 

The following theorem is a generalization of Theorem [l] to non- 
differentiable cases. Derivatives are replaced by generalized derivatives 
in a Chebyshev scale, and box dimension by similar notion of critical 
Minkowski order with respect to a Chebyshev scale. It shows that, 
in non-differentiable cases, the length of the e-neighborhood should be 
compared with the inverted Chebyshev scale instead of the power scale 
to obtain multiplicity. 

Let T = {f\\\ G A} be a family of functions on [0, d) admitting a 



uniform asymptotic development (11) in a family of Chebyshev scales 



T\ = (u (x, A), Ut(x, A), . . .), Q = (g x ) = {id - f x \fx G J 7 }. Let, for 
A = Ao, the monomials in the scale X = Xx be positive and strictly 
increasing on (0, d), for i > 1. 

Theorem 3. Let f = fx be a function from the family J 7 above, 
satisfying all assumptions of Theorem\^and the upper power condition. 
Let g = id — f . Then the following claims are equivalent: 

(i) A(/)(0) = 0, for i = 0, . . . , k - I, and D k {f){0) > 0, for some 
k > 1 (that is, f ~ Uk, for some k > 1), 



(ii) \A e (S9(x )\ 

(iii) m(g,T) = k. 



-i, 



'■k 



If moreover A C M, N , k < N, and the matrix (^ L )i=o...k-i, j=i...fc(Ao) is 
of maximal rank (i.e. equal to k), then (1), (2) or (3) is also equivalent 
to 

O) Vo x (g,G) = k. 

Without this regularity assumption, (i), (ii) or (iii) implies 

^l lx (g,Q) < k. 

On the importance of the upper power condition in Theorem [3j see 
Remark [5] in Appendix. 
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In the differentiable case, we notice that differentiation diminishes 
critical Minkowski order by 1. Let / G Diff r [0,d) be differentiable 
enough and suppose fio(f') > 1. Put g —id—f and h —id—f, then by 
Theorem [l] and Remark |2}ii) we have 

m(h,X) = m(g,I) — 1, 

where X = {1, x, x 2 , . . . , x r } is a differentiable Chebyshev scale. 

The same property is valid in the non-differentiable case when / 
has asymptotic development of some order in a Chebyshev scale, if 
the derivative is substituted by the generalized derivative in that scale. 
The following corollary is a direct consequence of Theorem [3] 

Corollary 1 (behavior of the critical Minkowski order under differen- 
tiation). Let X = {u ,Ui, . . . ,Uk} be a Chebyshev scale and let -Di(X) 
denote the Chebyshev scale of first generalized derivatives of I, that 
is, -Di(X) = {Di(ui), Di(u 2 ), ■ ■ ■ , Dx(uk)}- Suppose that f has an as- 
ymptotic development in the scale X of order k on (0, d) and let f and 
Di(f) satisfy assumptions from Theorem^ and the upper power con- 
dition. Let g =id—f, h =id—D\(f). It holds 

m(h,D 1 (Z))=m(g,T)-l. 

Finally, let us explain the gain of introducing the notion of critical 
Minkowski order with respect to some Chebyshev scale over the stan- 
dard box dimension. As we have mentioned before, the role of the box 
dimension and Minkowski content was to measure density of the orbit 
S 9 (xq) around the fixed point by determining the rate of growth of 
\A e (S 9 (xo)) \ as £ — > 0. If we consider, for example, functions fi(x) = x k 
and f2{x) = x h (— logx), k > 1, and compute standard box dimension 
of the orbits generated by g\ = id — fi and g 2 = id — f 2 , we get in both 
cases ^ 

dim B (g 1 ) = dim^^) = 1 - r- 

Thus box dimension is equal for two functions obviously distinct in 
growth. This is unnatural, the difference can be noticed in the upper 
(lower) (1 — 1/ A;)— Minkowski content A^ 1-1 ^, which is zero for g 2 and 
greater than zero for gx, thus signaling that the orbit generated by g\ 
has bigger density around than the one generated by g 2 . The reason 
lies in the fact that the box dimension and the Minkowski content are 
defined in a way that compares functions to power functions, but the 
functions f 2 (x) = x k (— log a?) are not visible in the scale of power func- 
tions. Precisely, by Theorem^J ^(^(xo))! — f 2 1 (^) does not behave 
as power of s, but satisfies e 1 ^ < \A e (S 92 (x ))\ < s l ^ k+& \ for every 5 > 
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0, as e tends to 0. Therefore there is not much sense in comparing the 
lenght to powers of e as in the standard definition of box dimension. We 
should define some other Chebyshev scale in which fi and f% both have 
developments, for example X = {1, x(— log a?), x, x 2 (— logx), x 2 , . . .}, 
and consider critical Minkowski orders with respect to this new scale 
instead of box dimensions. Then, we obtain distinct numbers: 



m(gi,I) = 2k, m(g2,Z) = 2k — 1. 

Therefore critical Minkowski order with respect to the appropriate scale 
is a more precise measure for density of the orbit in non-differentiable 
case then the box-dimension. 

3. Proof of the main theorems 

In the proof of Theorem [2] and of Theorem [3] we need the following 
lemma: 

Lemma 2 (Inverse property). Let d > and let f,g G C 1 (0,(i) be 
positive, strictly increasing functions on (0,d). 

i) If there exists a positive constant M > such that the upper 
power condition holds, 

(19) x ■ (log f)'(x) < M, x £ (0, d), 
then 

(20) (a) f^iy) — g 1 ^), as y — > implies f(x) ~ g(x), as x — > 0. 

fix) f~ x (v) 

(21) (b) lim if-r = (+oo) implies lim - i = +°° (°)- 

x^og(x) y^g^iv) 

ii) If there exists a positive constant m > such that the lower 
power condition holds, 

(22) m < x ■ (log f)'(x), x e (0, d), 
then 

(23) f(x)~g(x), as x — )■ 0, implies f~ l {y) — g~ (y), as y — > 0. 
Proof. 

i)a) From f^ 1 ~ g^ 1 we have that there exist constants A < 1, 
B > 1 and 5 > such that 

Ag-\y) < f-\y) < Bg-\y), ye (0,5). 

Putting x = g~ l (y) and applying / (strictly increasing) on the above 
inequality we get that there exists 5\ > such that 

(24) f(Ax)<g(x)<f(Bx), x G (OA). 
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For each constant C > 1 we have, for small enough x, 

log f(Cx) -log f(x) = (log/)'(0(C-l)x 
(25) < (log/)'(0(C-l)£, £G(z,Ca;). 



Combining (19) and (25), we get that there exist constants mp > 1 
and dc > such that 



(26) 



f(Cx) 

m 



< m c , x G (0, d, 



c) 



Now using property (26) and inequality (24), for small enough x we 
obtain 

i 

-f(x) < g{x) < m B f(x), 



my a 

i.e. f(x) ~ g(x), as a; — > 0. 



z) 6) Suppose lim x ^ ^fj = +oo- We prove that lim^o ^r^j = by 
proving that limit superior and limit inferior are equal to zero. Suppose 
the contrary, that is, 

/-%) 

liminf — — — = M, for some M > or M = oo. 
v-*° 9 \V) 

By definition of limit inferior, there exists a sequence y n — > 0, as n — > 
oo, such that 



(27) 



g~ l {yn) 



— > M, as n — >• oo. 



From (27) it follows that there exist no G N and C > such that 

(28) g~\y n )<cr\y n ), n>n . 

Now, as in i)(a), by a change of variables x n = g^ 1 (y n ), x n — > 0, and 
applying / (strictly increasing) on (28), we get 

m c g{x n ) > f{x n ), n > no, x n -)■ 0, for m c > 0, 

which is obviously a contradiction with linx^o ^j|y = +oo. Therefore 

liminf ^4 = 0. 

y^o g 1 (y) 

It can be proven in the same way that limit superior is equal to zero. 
Now suppose lim x _3.o = 0. In the same way as above, we prove 

^ 9{ x ) 

that lim^o f=r^j = 0. 

ii) It is easy to see by change of variables x = f~ l {y) that property 



(22 ) of / is equivalent to property ( 19 ) of / 1 and the statement follows 
from i). □ 
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Remark 3 (Counterexamples in Lemma p2J). To show that in Lemma 



ty i) the upper power condition (19) is important, we can take, for 



example, functions f(x) = e ^ and g(x) = e * . They do not satisfy 



(19) and, obviously, 



lin i4-T = 00 ' hut f 1 (y) = ~wT- — - 9 1 (y) - 



>o g(x) 2 logy ■ logy 



We can do the same for the lower power condition ( 22 ) in Lemma [^| 
ii) by considering, for example, f(x) = — and g(x) = — 2Ic ) gJ . - 

Proof of Theorem 

From lower power condition together with /'(0) = 0, we get that 
f(x) = o(x) and that f(x) is strictly increasing on (0, d). It can easily 
be checked that x n — > and d(x n ,x n+ i) — > 0, as n — > oo. Denote 
by N e and T £ nucleus and tail of the e— neighborhood of the sequence, 
that are e— neighborhoods of two subsets of the orbit satisfying the 
inequality d(x n ,x n+ i) < 2e for the nucleus, and d(x n , x n+ \) > 2e for 
the tail. Therefore, 

(29) \MS°( Xl )) = \N E (S 9 (x 1 ))\ + \T e (S°( Xl ))\, 

where |iV e | is the length of the nucleus, and \T e \ the length of the tail 
of the 5-neighborhood. For more on notions of the tail and the nucleus 
of the e— neighborhood of a set, see Tricot [16]. 

To compute the length, we have to find the index n £ G N such that 

(30) f(x ne ) < 2e, /(x ne _x) > 2e, 

that is, the smallest index n e such that e— neighborhoods of the points 
x Ue , x n e +ii e ^ c start to overlap. 
Then we have 

(31) |iV £ | = x ne +e, 

(32) \T £ \ ~ n £ -e. 

First we estimate \N £ \. From f(x) = o(x) we get 

(33) i im =Q 

w-+o/ 1 {y) 



Since / is strictly increasing, from (30) we easily get x Ue ~ f (2e). 
Since / satisfies the lower power condition, by Lemma |2}m) it follows 
x nE ~ / _1 (^)- This, together with (31) and (33), implies \N £ \ ~ f^ 1 ^)- 
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1/f(V> 



1/f(Xn D+ l) 




x n E [t 



Xr fo+1 d 
x np-t2 | 

AX n() + 1 AX,^ 



1/f(Xn c ) 



1/f ( x n +l) 

1/f(x no ) 




FIGURE 1. Sums from (1361) as sums of areas of rectangles. 



Now let us estimate the length of the tail, \T e \, by estimating n £ . 
Putting Ax n := x n - x n+1 , from x n+1 - x n = -f{x n ) we get 



(34) ^ = i and jr A ' ' 



f(Xr 



n=no 



1 = n £ — n ~ n £ , as e — > 0, 



n=no 



for some fixed n G N. 



As in (41) below, we get that tends to 1, as n tends to infinity, 



and thus we can choose the integer no so that 

(35) Af(x n+ i) < f(x n ) < Bf(x n+1 ), n > n , 

for some constants A, B > 0. 

Since the function jj^ is strictly decreasing on (0, d) and lim x ^ jk* 



+oo, the sum „ is equal to the sum of the areas of the rect- 

angles in Figure 1 i) and, analogously, the sum f , Axn ; is equal 

to the sum of the areas of the rectangles in Figure fljii). Therefore we 
have the following inequality: 



(36) 



Ax ra < T"o _cfe_ y> Ax n 

n=n J V ^We+l J V ^ n=nQ ^ V ""t-J-> 



E 
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From (35), we get 
(37) 



n=no 



AXr 



f(Xn 



< 



+1, 



n=no 



/Or 



so finally, putting (37) in (36) and using (|34j), we get the following 
estimate for n E : 

Xn ° dx 



(3 



, +1 /(*) 



as £ — > 0. 



Substituting x = f 1 {y), from the lower power condition we get 

(T7(y) 

— > m 

y 2 y 

and, consequently, for y G (0, f(d)), 

(39) _ f r%) y = _(T7M + /"'(v) 



y 



> (m - 1) 



Now substitution x = f 1 (s) in the integral (38) together with (39) 
gives 



(40) n e 
It holds 



CT7(g)jg < i 



m — 1 



/(*no) 
/(*n e +l) 



/fcnj _ /(«n e -l - /(»n,-l)) 



/On e -l) /On e -l) 

/(sn e -i) + /'&)(-/(sn.-i)) 

/(^n s -l) 

for some £ e G (x ne , x„ e _i), so /'(0) = implies 



1 -/'(&), 



(41) 



Hm 



1. 



Tie— l j 



From (30) and (41), we now conclude that f(x ne+ i) 
(40) becomes 



e. Therefore 



for some C > 0. From (32), we have that 

for some C\ > and e small enough. Together with |iV e | ~ f~ l {^) 
obtained above, this implies, using (29), that 

|^(5f(x 1 ))|~/- 1 (e), ase^O. 
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□ 



Proof of Theorem [3J 

We first prove that (i) =>• (ii) =>■ (Hi). Suppose Z)j(/)(0) = 0, 
i — 0, . . . , k — 1, .Dfc(/)(0) > 0, i.e., / ~ as x — >• 0. Theorem [^ap- 
plied to / gives \A e (S 9 (x\))\ — f~ x (e). Since / ~ Uk, by Lemma^n) 
we get f^ 1 ~ u^ 1 and therefore \A £ (S 9 (xi))\ ~ w^ 1 ^). Since Uk sat- 
isfies upper power condition, by Lemma [2] and Definition [5j of critical 
Minkowski order (see Remark j2jii)), we get m(g,X) = k. 

Now we prove that (Hi) =>• (ii) =>- (i). Suppose m(g,X) = k and 
/ ~ ui, for some I ^ k. As above, we conclude m(g,I) = I ^ k, which is 
a contradiction. Therefore f — Uk and, again as above, \A e (S 9 (xi))\ ~ 

By Lemmabl we conclude that (i) implies ^^(g, Q) < k. If moreover 
the condition on the maximal rank of the matrix (|f i )i=o...fc-i, j=i...&(Ao) 
is verified, by Lemma [l] («) is equivalent to (zu). 

□ 



4. Applications 

4.1. Cyclicity of limit periodic sets for planar systems. 

The number of limit cycles that bifurcate from a monodromic limit 
periodic set in an unfolding is equal to the multiplicity of the isolated 
fixed point x = of the Poincare map in the family of Poincare maps for 
the given unfolding, see e.g. Proposition 2 in |4]. For exact definitions 
of limit periodic set and cyclicity, see e.g. Roussarie [15]. 

Results from Section [2] connect cyclicity of a limit periodic set in an 
unfolding and the rate of growth of the length of the e-neighborhood 
of any orbit of the Poincare map whose initial point is sufficiently close 
to the limit periodic set. This rate of growth is given by the critical 
Minkowski order with respect to the appropriate scale. 

The orbit of the Poincare map on a transversal to the limit periodic 
set is the intersection of the corresponding one-dimensional orbit of the 
vector field with the transversal. Locally in a neighborhood of a point 
on the transversal, the structure of the one- dimensional orbit is that 
of the orbit of the Poincare map by a segment. Hence all interesting 
data of the one-dimensional orbit are given by the corresponding zero- 
dimensional orbit of the Poincare map. 

Therefore, instead of considering the rate of growth of the area (2- 
Lebesgue measure) of the e-neighborhood of an orbit of the field itself, 
as e — > 0, which would be more natural in search of cyclicity of a limit 
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periodic set, it is sufficient to consider the rate of growth of the length 
of the e-neighborhood of an orbit of its Poincare map. 

Let T be the stable limit periodic set for the analytic unfolding X\. 
We consider only the unfoldings of finite codimension such that the 
family of Poincare maps gx(x), x G [0,d), for the unfolding is well 
defined and different from identity on the transversal to vector field X\ 
in neighborhood of V. 

Let us recall that the function f\=id—g\ is called the displacement 
function. The main idea is to find the family of Chebyshev scales T\ 
such that the family f\ has a uniform development of some order in the 
family Z\, as was introduced in Section [2] Suppose T is a limit periodic 
set of X Xo , for the parameter value A = A . Then, by Theorem |3j the 
critical Minkowski order of g\ Q with respect to the scale Z\ , m(gx , I\ ), 
is an upper bound on cyclicity of T in the unfolding X\. 

In the sequel, limit periodic sets are stable limit cycles, nondegen- 
erate stable focus points and stable homoclinic loops. The family of 
displacement functions fx for the unfolding has a uniform asymptotic 
development in a family of appropriate Chebyshev scales and is analytic 
in x = in first two cases ( different iable cases) and non-differentiable 
in x = in the case of homoclinic loop. 

4.1.1. Differentiable case, limit cycle. 

Suppose that X\ has a stable or semistable limit cycle T and let X\ 
be an arbitrary analytic unfolding of X\ . 

There exists neighborhood W of Ao such that the displacement function 
fx is analytic on [0, d), for A G W, and /a q (0) = 0. Expanding f\(x) in 
Taylor series, we get 

(42) f x (x) = a (\) + a x (X)x + a 2 (X)x 2 + a 3 (X)x 3 + . . . , A G W. 

The family f\ has a uniform asymptotic development in the Cheby- 
shev scale 

X {1, x, a; , . . . , x }, 

of any order I G N. 

By Theorem [2j the length of the e-neighborhood of an orbit of the 
Poincare map around the limit cycle should be compared to the inverted 
scale, {e, e 1//2 , e 1 ^ 3 , . . .}, to obtain an upper bound on the cyclicity. Let 
Q = {g\ = id — fx} be the family of Poincare maps. By Theorem [3j 
if fx ~ x k , for some 2 < k < £, then \A S \ ~ e 1//fe and the critical 
Minkowski order of g\ with respect to T\ is equal to k, m(gx ,T\ ) = 
k. The cyclicity of the limit cycle in the unfolding Xx is equal to 
{9\oiG) — k- ^ moreover the unfolding (Xx) is general enough 
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so that the regularity condition from Theorem [3] is satisfied, then the 
cyclicity /j, f lx (g Xo ,G) = k. 

4.1.2. Differentiable case, weak focus. 

Suppose Xq IS £L stable weak focus of X Ao (that is, DXx { x o) has two 
conjugate complex eigenvalues without the real part). Suppose X\ is 
an arbitrary analytic unfolding of X\ . 

There exists neighborhood W of Ao such that, for A G W, the dis- 
placement function f\(x) is analytic in and f\(0) = 0. Therefore we 
can expand fx in Taylor series around and, by symmetry argument 
around focus point, we get that the leading monomials can only be the 
ones with odd exponents: 
(43) 

f x (x) = (3 1 {X){x+g 1 {X,x))+(3 3 {X){x i +g 3 {X,x))+/3 5 {X){x 5 +g 5 {X,x))+. . . , 

where gi(X, x) denotes some linear combination of monomials from Tay- 
lor expansion of order strictly greater than x l and with coefficients 
depending on A. 

The family of displacement functions fx has a uniform asymptotic 
development in a family of Chebyshev scales Ix of some order: 

1\ = {x + gt(X,x), x 3 + g 3 (X,x), x 5 + g 5 (X,x), . . . , x 2e+1 + g 2 £+i(X,x)}. 

To obtain an upper bound on the cyclicity of the focus, by Theo- 
rem [3j the length of the e-neighborhood of the discrete orbit of the 
Poincare map around the origin should be compared to the inverted 
scale {e, e 1 ^ 3 , e 1 ^ 5 , . . .} of Z\ . We proceed as in the example above. 

4.1.3. Non- differentiable case, homoclinic loop. 

Suppose Xx has a stable homoclinic loop with the hyperbolic saddle 
point at the origin, with Ox as unstable and Oy as stable manifold, and 
such that the ratio of hyperbolicity of the saddle is r(A ) = 1 (i.e. 
_DX Ao (0) has eigenvalues of the same absolute value, but of different 
sign). Suppose X x is an analytic unfolding of X Xq and that, for A G W, 
each Xx has a hyperbolic saddle of ratio r(A) at the origin, with the 
same stable and unstable manifolds. 

We consider the family of Poincare maps Q = g\ and the family of 
displacement functions fx = id— g\, x G (0, 5), on a transversal to stable 
manifold near the origin, as in Chapter 5 in Roussarie [15J. The family 
cannot be extended analytically to x = due to nondifferentiability in 
x = and the following asymptotic expansion in x = holds instead 
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(see [15], Section 5.2.2): 

fx(x) = /3 (\) +a 1 (\)[xuj(x,ai(\)) + gi{x,\)} + 

+ Pi{X)x + a 2 {X)[x 2 u{x, ai(A)) + g 2 (x, A)] + (3 2 {\)x 2 + ...+ 
+ /3 n (A)x + a n (A)[x n w(x,a 1 (A)) +g n (x,X)] + f3 n (\)x n + o(x n ), 
(44) nGN, 

where an (A) = 1 — r(A), gi(x, A) denotes linear combination in the 
monomials of the type x k oj l of strictly greater order than x l u (order on 
monomials is defined by increasing flatness, x l uji < x k u l if (i < k) or 
(i = k and j > I)) and 



u(x, a) 



if a ^ 0, 



— logx if a = 0. 

The family of displacement functions f\ has obviously an uniform 
asymptotic development in the following family of Chebyshev scales T\ 
of some order: 

Z A = {l,xu(x, oii(X)) + gx{x, X),x,x 2 u(x, «i(A)) + g 2 (x, A), x 2 , . . .}. 
If we take A = A in (44), we get the following expansion for X\ 



(ai(Ao) = 0, /a (0) = 0): 

fx (x) = 0i{X o )x + a 2 (\ )x 2 u(x,0) + a 3 (\ )x 3 u(x,0) + ... = 

(45) = /3i(A )x + a 2 (A )x 2 (-logx) + a 3 (A )x 3 (- logx) + ... 

The length of the e-neighborhood should be compared to the inverted 
scale of I\ to obtain information on cyclicity. The critical Minkowski 
order signals the moment the comparability occurs. By Theorem |3j if 
fx (x) ~ x k , as x — > 0, k > 2, then the critical Minkowski order is 
equal to 2k, m(g\ ,I\ ) = 2k; if f\ ~ x k (— logx), k > 2, then the 
critical order is equal to 2k — 1, m(g\ ,T\ ) = 2k — 1. Consequently, 
the cyclicity of the loop is less than or equal to 2k, 2k — 1 respectively. 
Equality can be obtained if the unfolding (X\) is general enough so 
that the regularity condition from Theorem [3] is satisfied. 

4.1.4. Non-differentiable case, Hamiltonian hyperbolic 2-cycle with con- 
stant hyperbolicity ratios. 

Suppose (X\) is an unfolding of a Hamiltonian hyperbolic 2-cycle of 
the field X\ in which at least one separatrix remains unbroken. Such a 
situation appears for polycycles having part of the line at infinity as the 
unbroken separatrix. Suppose that the ratios of hyperbolicity of both 
saddles Si and S 2 at A = A are ri(Ao) = ^(Ao) = 1- The breaking 
parameter of the breaking separatrix is denoted by /?i(A) (/?i(A ) = 0). 
By x G (0, 5) we parametrize the (inner side) of the transversal to the 
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stable manifold of one of the saddles, and we choose the saddle whose 
stable manifold is on the unbroken separatrix, say Si. In search of 
cyclicity, instead of considering fixed points of Poincare maps f\ on 
(0, S), for simplicity we can consider zero points on (0, 5) of the family 
of maps 

A A (x) = D^o R^(x) -R{o D$(x), 

where D 1 and D 2 represent Dulac maps of the saddle Si, R\ is the reg- 
ular map along the broken separatrix and R 2 the regular map along the 
unbroken separatrix. Obviously, -Ri(O) equals the breaking parameter 
of the separatrix, /?i(A), and R 2 (0) = 0, on the unbroken separatrix. 
Using the developments of Dulac maps from [15], just like in the above 
example of the saddle loop, Aa has the uniform development in the 
monomials from the two Chebyshev scales, X\ and X\ below, since the 
developments for D\ o R\{x) and R\ o D^x) are subtracted: 

X\ = {1, xui(x, «i(A)), x, x 2 u\ (x, «i(A)), x 2 oji(x, ai(A))), x 2 , 
x 3 u>l(x, ai(A)), x 3 tu 2 (x, «i(A)), x 3 wi(x, «i(A)), x 3 , . . .}, 

X\ = {l,xu 2 (x,a 2 (X)),x, x 2 uj\{x 1 a 2 (X)), x 2 u 2 (x, a 2 (A))),x 2 , 
x 3 cu 2 (x, a 2 (X)),x 3 u 2 (x, a 2 (X)),x 3 u 2 (x, a 2 (X)),x 3 , . . .}. 

For the development, see e.g. [2]. For each monomial x k u\ 1 k > 1, / > 
0, it necessarily holds that k > I, «i(A) = 1 — ri(A), a 2 (X) = 1 — r 2 (A), 
and cui and u 2 are as defined in the section above. They are known as 
independent compensators, since they are not comparable by flatness, 
and thus disable the concatenation of X\ and X\ in one Chebyshev 
scale. 

If we additionally suppose that the ratios of hyperbolicity ri(X ) and 
r2(Ao) are preserved throughout the unfolding, then we have 

u>i(x, ai(A)) = u 2 (x, a 2 (A)) = — logx, for all A. 

In this case the Chebyshev scale in which all of from the unfolding 
{X\) have the uniform development is 

X ={1, x, x 2 {— logx) 2 , x 2 {— logx), x 2 , 

x 3 (— logx) 3 , x 3 (— logx) 2 , x 3 (— logx), x 3 , . . .}. 

To see the number of limit cycles that can arise in the unfolding of the 
hyperbolic 2-cycle in X\ , by Theorem [3j the length of e-neighborhood 
of the discrete orbit of A\ (x) should be computed numerically and 
compared to the inverted scale of X. The index i for which \A £ \ ~ 
u^ 1 {e) holds, in the article called critical Minkowski order m(A Ao ,X), 
represents an upper bound on the number of limit cycles that can 
appear in the unfolding (X A ) of X\ . 
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Let us note here that this upper bound is not necessarily optimal, 
since the scale X is taken to be the largest possible for a given problem. 
Better results on upper bound are obtained in [I], using asymptotic 
developments of Abelian integrals, and in [8]. In [8], the upper bound 
is given in terms of characteristic numbers of holonomy maps, not using 
asymptotic development of the Poincare map. 

4.2. Abelian integrals. 

Abelian integrals on 1-cycles are integrals of polynomial 1— form u 
along the continuous family of cycles of the polynomial Hamiltonian 
field, lying in the level sets of the Hamiltonian H, 5t C {H = t}, 



Suppose that the value t = is a critical value for the Hamiltonian 
field in M 2 , such that there exists d > and a continuous family of 
cycles belonging to the level sets {H — t}, t G (0,<i). Then we have 
the following asymptotic expansion at t — (see Arnold pQ, Ch. 10, 
Theorem 3.12 and Zoladek [T7], Ch. 5): 



where a runs over an increasing sequence of nonnegative rational num- 
bers depending only on Hamiltonian H (x, y) (such that e 2ma are eigen- 
values of monodromy operator of the singular value) and a^^ G M.. 

Obviously, the corresponding Chebyshev scale for this problem is: 



(47) X = {t ai ( - log t) , t ai , t a2 ( - log t ) , t a2 , . . . , t am ( - log t) , t am , . . . } . 



It makes sense to compute critical Minkowski order of the orbit 
S 9 (xo), comparing the length of e-neighborhood of g(t) = t — I Uo (t) 
with the inverted scale of X, to obtain the multiplicity of an Abelian 
integral in a family of integrals. 

In IR 2 , Abelian integrals have been used as a tool for determining 
cyclicity of vector fields, considering them as perturbation of Hamil- 
tonian field (for details and examples see e.g. Zoladek [T7] . Ch. 6). 

Suppose we have the following rj— perturbed Hamiltonian system, 




i 



(46) 




a k=0 



(48) 



x 



y 



dH 

dy 
dH 

dx 



+ r)P(x,y,r)), 



+ r]Q(x,y,r]) 
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where P, Q, H are polynomials and r\ > 0. 

Let u = Qdx — Pdy be the polynomial 1— form defined by P, Q. 

Let t = be a critical value of Hamiltonian, and let S be a transversal 
to the family of cycles (St C {H = t}) on small neighboourhood of 
t = 0, parametrized by t G [0,d). Then (see e.g. Zoladek [17], Ch. 6) 
the displacement function on S of the perturbed Hamiltonian field is 
given by 

(49) f v (t) = V I u (t) + o( V ), 

i.e., Abelian integral is the first approximation of the displacement 
function, for 7] small enough. Here we suppose that I w (t) is not identi- 
cally equal to zero, i.e. that u is not relatively exact. 

On some segment [a, (3} C (0, d) away from critical value t — 0, it 
is known that the number of zeros of Abelian integral gives an upper 
bound on the number of zeros of the displacement function f v (t) on 



a,/3] of the perturbed system (48), for rj small enough (both counted 



with multiplicities), i.e. on the number of limit cycles born in perturbed 
system (48) in the area Ute[a/3]^*' ^ or V < Vo small enough (for this 



result, see e.g. [3], Theorem 2.1.4). 

However, the problem arises if we approach the critical value t = 
and the result cannot be applied to the whole interval [0,d). In some 
systems, some limit cycles visible as zeros of displacement function are 
not visible as zeros of corresponding Abelian integral, because some- 



times the approximation (49) is not good enough. One of the examples 
is the perturbation of the Hamiltonian field in the neighborhood of the 
saddle polycycle with 2 or more vertices, see Dumortier, Roussarie jl]. 
Abelian integrals near hyperbolic polycycles have an expansion linear 



in logt, see expansion (46) or [4J, Proposition 1. On the other hand, 
see Roussarie [15] , the asymptotic expansion of the displacement func- 
tion near the saddle polycyle with more than one vertex involves also 
powers of logt greater than 1. 

In the neighborhood of the center singular point and of the sad- 
dle loop (1— saddle polycycle) of the Hamiltonian field, however, the 
multiplicity of corresponding Abelian integral gives correct information 
about cyclicity, see e.g. Dumortier, Roussarie [4], Theorem 4. 



4.3. One example out of scope of Theorem [3[ 

At the very end, let us note that in the former examples we have 
used critical Minkowski order which reveals the rate of growth of e- 
neighborhood of the orbit generated by Poincare map around limit 
periodic set to conclude about the cyclicity of the set. The connection 
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is given by Theorem [3] through the notion of multiplicity of the fixed 
point zero which is equal to cyclicity. 

From the assumptions of Theorem [3] it is visible that the theorem 
cannot be applied to displacement functions which are infinitely flat 
and therefore not comparable to powers (see the paragraph after Defi- 
nition [4]). We have noticed that this restriction of Theorem [3] to func- 
tions that are not infinitely flat makes sense in applications. As an 
example, we can take the accumulation of limit cycles on the saddle- 
node polycycle, a case which obviously should not meet the conditions 
of Theorem [3] because multiplicity and cyclicity should not turn out 
finite. It is interesting that in this case the displacement function is 
infinitely flat, and therefore excluded from Theorem [3j Perhaps this 
could be the subject of further research. 

We state here the result from Il'yashenko [10J: If a sequence of 
limit cycles of an analytic vector field converges to a polycycle with 
saddle-node singular points, then one can select a semitransversal to 
this polycycle such that the displacement function is not equal to zero, 
but infinitely flat, for e.g. f(x) ~ e~*. If we have a polycycle with 
only saddle singular points, then the displacement function cannot be 
infinitely flat. 



5. Appendix 

In Appendix we put some observations concerning main results. 



Remark 4 (sublinearity in Theorem |2J). The condition m > 1 in the 
lower power condition in Theorem [| cannot be weakened. If we take, 
for example, the function 

m 



— log X ' 



it obviously satisfies all assumptions of Theorem [§| except sublinear- 
ity: the lower power condition holds only for m < 1. If we compute 
\A e (S B (xo))\ for this function (as it is computed in the proof of Theo- 



rem 



in Section 



Sp ; we get that l^j^gHJ tends to infinity, as e — > 0, 



conclusion (17) is not true 



and therefore the 

On the other hand, for functions of the form 

f(x) = — , a > 0, 

— log a; 

which are obviously sublinear with m = 1 + a > 1, the explicit compu- 
tation shows that \ A e (S 9 (xq))\ — f^ 1 ^), as e — > 0. 
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Remark 5 (Upper power condition in Theorem [3]) . The upper power 
condition on f is needed in Theorem^ for the if implication to hold, see 
Lemma^i). As a counterexample, we can take the following Chebyshev 
scale 

X = {e", , e~a^ , . . .}, 

and, f(x) = e~^, g=id—f , which does not satisfy upper power condi- 
tion. 

Obviously, D (/)(0) = £>i(/)(0) = and D 2 (f)(0) > 0, so the mul- 
tiplicity Ho(g,G) < 2. On the other side, (e) ~ (e) — u^ 1 (e) ~ 
... ~ _ log e > therefore the critical Minkowski order m(g,Q) is infi- 
nite. In this case, we are not able to read the multiplicity neither 
from the critical Minkowski order nor from behavior of the length of 
the e— neighborhood. 



Example 3 (Non-flat, non weakly comparable to powers function). 
We construct a non infinitely flat function f that does not satisfy 
x{\ogf)'(x) < M for any M > 0, just to show that, for functions of 
interest, non-flatness is not equivalent to weak comparability to powers. 

The main idea is to bound the function from above and from below 
with x a+1 and x a , a > 1, therefore it cannot be infinitely flat. 
Next we need to make sure that on some intervals approaching zero its 
logarithmic growth is faster than the logarithmic growth of x a . 

We construct the function f in logarithmic chart, i. e. we construct 
function g(x) = \ogf(x) on some segment (0, d). 
Let hi(x) = log(a; a ) = a logs and let h 2 (x) = log(a; a+1 ) = (a+1) \ogx. 
Let us take X\ close to x = 0. The segment Ii connects the points 
(xi,hi(xi)) and (xi/2, h 2 (xi/2)) . Now we choose point x 2 such that 
hi(x 2 ) < h 2 (xi/2) (to ensure that f is increasing). We get segment I 2 
by connecting (xi/2,h 2 (xi/2)) and (x 2 ,hi(x 2 )). We repeat the proce- 
dure with x 2 instead of x\ to get segments 1% and 1^ and, inductively, 
we get the sequence (x n ) tending to 0, as n — > oo, and the sequence 
of segments (I n ) which are becoming perpendicular very quickly, see 
Figure \^ 

The graph of our function g will be the union of the segments U^Li In, 
smoothened on edges. Obviously f(x) = e 9 ^ is bounded by x a+l and 
x a . Nevertheless, if we take the sequence (y n ) such that x n /2 < y n < 
x n , we compute 

mlogx n - (m + 1) log ^ 
9 [yn) ■ Vn = -py Vn — - logx„, as n ->• oo, 
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x,/2 




FIGURE 2. Function g(x) = logf(x) from Construction [3] 

and, thus, for the sequence (y n ) tending to it holds g'{y n )yn ~^ oo? 
n — > 00, a contradiction to xg'(x) < M. 

Acknowledgements: We would like to thank Darko Zubrinic for 
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